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Abstract 

We investigate island formation during heteroepitaxial growth using an atom- 
istic model that incorporates deposition, activated diffusion and stress relax- 
ation. For high misfit the system naturally evolves into a state characterized 
by a narrow island size distribution. The simulations indicate the existence 
of a strain assisted kinetic mechanism responsible for the self-assembling pro- 
cess, involving enhanced detachment of atoms from the edge of large islands 
and biased adatom diffusion. 
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Heteroepitaxial growth of highly strained structures has gained interest lately as it offers 
the possibility to fabricate islands with very narrow size distribution without any substrate 
patterning process. Thanks to their small size, these islands, coined self-assembling quantum 
dots (SAQD), are candidates for one dimensional electron confinement SAQD formation 
has been observed for a wide range of material/substrate combinations, including InAs on 
GaAs U-U , InGaAs on GaAs 0J§, AlInAs on GaAlAs §, GeSi on Si HID], InP on InGaP 



121 , GaSb on GaAs [13| and ZnSe on ZnMnSe [|T4 | , indicating the existence of a not yet 



understood common mechanism governing the self-assembling process. 

In this paper we investigate the kinetics of island formation during heteroepitaxial growth 



with lattice constant a° on a substrate with lattice constant a°, we find that sufficiently large 
misfit, e = (a * — a°)/a°, leads to self-assembled island formation in the system. In particular, 



using a one dimensional model [15j that includes all microscopic elements common to the 
materials for which SAQD formation has been observed, namely deposition, activated diffu- 
sion, and strain relaxation at every deposition and diffusion event. Depositing 2ML atoms 

° on a substrate with lattice constant a°, 

/-«*.; ^ 

e = 5% and 7.5% leads to a narrowly peaked island size distribution, in contrast with a wide 
distribution for e = 0% and 2.5%. We discuss the kinetic mechanism responsible for the self- 
assembling process, and compare our results with experimental work on SAQD formation. 
Since the model does not contain material dependent features, the mechanism responsible 
for self-organization is expected to be generic, applicable to a wide class of materials. 

Monte Carlo Method and Stress Relaxation — To include stress in the model, we con- 
sider that the atoms interact harmonically with their nearest and next nearest neighbors 
| 16| . The elastic strain energies are recalculated after every deposition or diffusion event, 
using a checkerboard relaxation method. The relaxation starts at the last active atom, and 
propagates radially outward, until in a two particle wide shell all relative displacements are 
smaller than a preset parameter 5. The relaxation is iteratively restarted from the same 
origin until even in the closest shell the displacement is smaller than 5. Radial relaxation is 
more efficient for our problem than algorithms relaxing the entire system, since most changes 
in the displacement occur where a particle moves, and decrease fast with the distance from 



the source. Trial runs with different values of 5 indicated that for 5 < 10~ 3 the results were 
practically identical. In the simulations we used 5 = 10~ 4 . 

The surface particles are allowed to hop to neighboring lattice sites, with an SOS condi- 
tion, disallowing up or down jumps larger than one atom high. The hopping probability for 
an atom is proportional to exp[— (nE n + E — E^j/ksT], where n is the number of neighbor- 
ing atoms, E n is the bond energy and Eq is the diffusion barrier for an isolated atom on a 
stress free substrate []l7f . The strain energy is given by E s = (c/2) XX a « — a ? )V ( a i > ) 2 where 



af and ai are the bulk and stretched bond lengths, i running over the occupied nearest 
and next-nearest neighbors. We use c = 44eV for the force constant, a typical value for 



many semiconductors fl8fl , E n = 0.3eV, Eq = 0.4eV and T=800K. The substrate consists 
of N = 50ML atoms with lattice constant a Q s = 1, on which we deposit with a constant 
flux 2ML atoms with lattice constant a® — a° s (l + e). The system size is L = 200 fT9| . We 
identify as island every mound with height larger than 1ML, and define the base size of the 
island, denoted by s, as the lateral size of the island measured in the second monolayer (to 
distinguish it from the wetting layer). The islands are coherently strained and dislocations 
are not allowed. 

Numerical results — The most convincing evidence of the stress induced self-assembling 
process is provided by the island size distribution, shown in Fig. |l|a. For e = 0% and 2.5% 
the distribution is wide, i.e. the system contains islands of all sizes, with a small peak around 
s = 20. However, for e = 5% and 7.5% the distribution has a narrow peak centered at s = 6 
for e = 5% and s = 5 for e = 7.5%. 

The parameter capturing the dynamics of self-assembly in the system is the relative width, 
w s /s, shown in Fig. 0a, where w 2 s = s 2 — s 2 is the width of the island size distribution and s 
is the average island size. An increasing w s /s indicates unbounded growth of fluctuations, 
while a decreasing one is a signal of self-organization in the system. As Fig. |I|a indicates, for 
e = 0% and 2.5% w s /s increases continuously with coverage, while for e = 5% and 7.5% w s /s 
increases only until it reaches a peak at some small coverage 8 C , after which it decays. The 
peak signals the onset of self-organization: for 9 > 6 C we witness a continuous increase in 



the uniformity of the island size. The peak is at 8 C = 0.66ML for e = 5%, and 9 C = 0.15ML 
for e = 7.5%, indicating that the self-assembling process is more effective for larger misfit. 

An experimentally often measured quantity is the island density, p. As Fig. [l|b indicates, 
for the stress free system p has a peak at 1.33ML, after which it decreases. This behavior is 
characteristic for homoepitaxy [F2TJ: for small coverages the dynamics is dominated by island 
nucleation. After a certain p is reached the incoming atoms are captured by the existing 
islands, prohibiting further island nucleation, and stabilizing the island density. Continuing 
the deposition leads to island coalescence, that results in a drop of the island density. Fig. 
[l]b is consistent with this scenario for e = 0% and e = 2.5%. However, we observe no 
such peak for e = 5% and 7.5%, indicating continuous island nucleation, without significant 
coalescence |2(| . 



Mechanism of self- organization — The main difference between the stress free and the 
stressed system comes in two strain related effects, that we discuss separately. 

(a) Strain lowers the energy barrier for diffusion, thus making diffusive hops more prob- 
able. Fig. |2] shows the strain energy in the vicinity of an island for e = 7.5%, indicating that 
the substrate is strained and that E s decreases as we move away from the edge of the island. 
This means that if atoms are deposited near the island, strain biases their otherwise random 
motion, generating a net surface current j = — V/i(x), where p(x) is the local chemical po- 



tential pl|| . The chemical potential is p ^ —(nE n + E — E s ), where nE n + E s is independent 
of the atom position for an isolated adatom on a flat surface. The only contribution to the 
current comes from the position dependence of the strain energy, leading to j ~ —VE S , 
that points towards the decreasing strain direction. Thus the strain field around an island 
generates a net current of adatoms away from the island. 

(b) For large islands the strain energy, E s , at the edge becomes comparable to the bonding 
energy of the edge atom, nE n + E (with n = 1), enhancing its detachment, thus leading 
to a gradual dissolution of the island. Such mechanism favors a smaller average island size 
and leads to a narrower island size distribution, as observed by Ratsch et al. [22|. 



The simultaneous action of (a) and (b) leads to a kinetic mechanism stabilizing the 



island size: as islands grow, a strain field develops, that helps to dissolute the edge atoms 
(effect (b)) and "pushes" them away from the islands (effect (a)). Furthermore, the newly 
deposited atoms also diffuse away from the larger islands (effect (a)). These combined effects 
slow the growth rate of large islands and increase the adatom density away from them, thus 
enhancing the nucleation of new islands. The newly nucleated islands are small, and so is 
the strain field around them, thus they grow at a much faster rate than the older and larger 
one. This eventually leads to a narrow island size distribution in the system (Fig. [l|a). The 
final island size is determined by a dynamical balance between the adatom density, the strain 
induced current away from the islands and the strain energy of the island, governing the 
detachment of the edge atoms. If we could monitor in real time the growth, we would witness 
a continuous nucleation of islands, such that small islands grow fast, and stop growing after 
they reach a certain size. 

Discussion — The observed behavior compares favorably with the main features of the 
experimentally observed SAQD formation. First, TEM observations of the strained islands 
document the existence of the strain field in the substrate . As Fig. |^ shows, such a 

field is reproduced by our relaxation method, and is responsible for the current j discussed in 
(a). Second, it is experimentally established that the increasing coverage contributes mainly 
to an increase in the island density, and much less to the further increase in the size of the 
existing islands @, which is reproduced by the simulations (see Fig. |l]b). Third, experiments 
on GalnAs growth on GaAs substrate document an increasing width in the early stages of 
the deposition process, followed by a gradual decay for larger coverages [@,|J. This is similar 
to the behavior shown in Fig. |l]a: w/s decreases only after a certain coverage 8 C has been 
reached. Indeed, for small coverages the islands are both small and distant, thus the strain 
induced biased diffusion (a) and atom detachment (b) are not yet relevant, and the island 
formation essentially follows a strain-free path. Only when 9 approaches 6 C the discussed 
strain-induced mechanisms reverse the growth in the relative width. Finally, the simulations 
indicate ordering in the distances between the islands, as observed for high island densities 
in some investigations |J. 
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I have benefited from enlightening discussions with J.K. Furdyna, M. Krishnamurthy, 
N.N. Ledentsov, J.L. Merz, M.S. Miller, K. Newman, and W. Seifert. 
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FIGURES 

FIG. 1. (a) Island size distribution measured after the deposition of 2ML atoms. Inset: 
Relative width w s /s as a function of coverage, (b) Island density as a function of coverage. In all 
figures the symbols correspond to misfit values 0% (o), 2.5% (□), 5% (O), and 7.5% (A). 

FIG. 2. The strain energy around a typical island. The substrate (filled □) and the islands on 
top of it (empty □) are shown on the upper part of the figure. E s is the strain energy of an atom 
placed on the top of the substrate or on the island. For example, E s at x = 18 is the strain energy 
felt by the adatom shown by the circle on the top of the substrate. One can see that E s is the 
largest when the atom is at the edge of the island (x = 21, 30). E s decays as the adatom moves 
away from the island, generating a net current, j(x), shown by the arrows. Note that E s does not 
decay to zero, since the monomer can locally stretch the substrate. 
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